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ABSTRACT 


Theoretical  methods  are  presented  for  the  analytic  and 
experimental  determination  of  mechanical  impedance  and 
its  associated  effects  on  the  dynamic  response  of  linear, 
multi-degree  of  freedom/multi- dimensional  systems. 
Impedance  matrices  are  formulated  which  are  useful  in 
describing  and  predicting  the  frequency  dependent  charac¬ 
teristic  properties  of  these  systems.  The  problem  of  inter¬ 
connecting  subsystems  is  formulated  in  terms  of  composite 
impedance  matrices  and  interconnection  equations,  from 
which  critical  frequencies  (resonances  and  anti- resonances) 
and  dynamic  response  of  the  composite  system  may  be  ob¬ 
tained.  Analytical  and  experimental  procedures  for  the 
evaluation  of  these  matrices  are  also  discussed.  To  demon¬ 
strate  the  application  of  some  of  the  procedures,  an  ex¬ 
ample  problem  is  solved.  Moreover,  the  results  of  this 
study  are  sufficiently  general  for  extending  the  approach  to 
many  types  of  distributed  system  configurations  subject  to 
varied  excitation  phenomena. 
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MECHANICAL  IMPEDANCE  ANALYSIS 
FOR  LUMPED  PARAMETER  MULTI-DEGREE 
OF  FREEDOM/MULTI-DIMENSIONAL  SYSTEMS 

by 

Frank  J.  On 

Goddard  Space  Flight  Center 


1.  INTRODUCTION 

The  analysis  of  mechanical  systems  under  dynamic  loadings  requires  that  the  dynamic  char¬ 
acteristics  of  the  systems  be  known  in  addition  to  the  characteristics  of  the  excitations.  In  com¬ 
plex  linear  systems,  it  is  useful  to  specify  these  system  dynamic  characteristics  in  terms  of 
mechanical  impedance  (References  1  and  2). 

In  this  study  a  mechanical  impedance*  analysis  of  complex  systems  is  developed.  This  anal¬ 
ysis  ultimately  may  serve  as  an  analytical  tool  in  the  computation  of  impedances  and  as  a  guideline 
in  the  laboratory  measurement  of  impedances. 

The  basic  philosophy  adopted  in  this  analysis  is  based  on  the  idea  of  breaking  up  a  complex 
system  into  component  parts  with  simple  dynamic  properties  of  mass,  stiffness,  and  damping  that 
can  be  formulated  readily  in  matrix  notation.  Appropriate  component  matrices  are  then  considered 
as  building  blocks,  that  when  fitted  together  according  to  a  set  of  predetermined  rules,  provide  the 
mechanical  impedance  characteristic  of  the  entire  system.  Such  a  philosophy  has  been  practiced 
in  the  well  known  method  of  influence  coefficients  (and  related  methods)  of  structurail  analysis.  The 
matrix  formulation  of  mechanical  impedance  may  be  considered  as  a  natural  extension  of  these 
methods,  and  like  them  it  is  well  suited  for  processing  by  digital  computers. 

The  results  of  these  extensions  to  the  determination  of  mechanical  impedance  require  the  re¬ 
interpretation  of  matrix  parameters  in  terminologies  that  can  be  associated  with  the  basic  concepts 
of  impedance.  This  is  accomplished  by  manipulating  equations  of  motions  to  appropriate  forms 
physically  corresponding  to  impedance  quantities  such  as  point  and  transfer  impedances,  free 
velocity  and  blocked  force,  etc. 


•Although  mechanical  impedance  and  mobility  are  respectively  defined  as  the  complex  ratio  of  sinusoidal  force  to  sinusoidal  velocity 
and  of  sinusoidal  velocity  to  sinusoidal  force,  four  other  ratios  between  sinusoidal  force  and  sinusoidal  displacement  and  acceleration 
are  equally  useful.  Xhus  the  developincnts  here  may  refer  to  all  six  ratios. 
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In  order  to  present  a  method  of  analysis  that  is  not  confined  to  any  particular  type  of  struc¬ 
tural  system,  a  general  method  based  on  multi-degree  of  freedom,  multi-dimensional  systems*  is 
adopted.  The  necessity  of  treating  real  complex  structural  systems  as  multi-dimensional  systems 
has  been  recognized  among  environmental  engineers  for  many  years,  but  impedance  concepts  for 
multi-dimensional  systems  have  been  advocated  only  recently.  This  is  largely  due  to  the  recogni¬ 
tion  that  as  aerospace  structural  systems  become  larger  and  more  complex,  one -dimensional  im¬ 
pedance  concepts  are  inadequate  (Reference  3).  Also,  as  the  demands  for  larger  environmental 
test  facilities  to  accommodate  these  systems  become  more  impractical,  the  use  of  subassembly 
and  component  testing  philosophies  becomes  imperative.  It  is  hoped  that  an  appreciation  of  the 
generalized  mechanical  impedance  approach  will  be  gained  from  the  developments  herein.  In  par¬ 
ticular,  two  aspects  of  the  approach  will  be  emphasized:  (a)  it  provides  both  a  means  of  analyzing 
complex  structural  systems  by  partitioning  into  successive  smaller  systems,  and  a  natural  means 
for  matching  boundary  conditions;  and  (b)  it  provides  a  useful,  sound  basis  for  performing  and  cor¬ 
relating  subassembly  and  component  environmental  tests,  obviating  the  requirement  for  huge  test 
facilities. 

In  the  following  sections,  the  theoretical  relations  and  methods  of  evaluating  the  pertinent 
parameters  required  for  the  successful  application  of  this  concept  are  presented. 


2.  EQUATIONS  OF  MOTION 


It  is  well  known  that  complex  linear  time  invariant  mechanical  systems  can  be  analyzed  in 
terms  of  their  lumped  parameter  equivalents.  The  matrix  equation  of  motion  for  such  a  grid  net¬ 
work  of  lumped  elements  is  given  by: 


NN  N1  NN  NI  NN  N1  N1 

where  x  denotes  generalized  coordinates  and  N  denotes  the  number  of  degrees  of  freedom.  The 
description  of  a  particular  system  is  contained  in  the  coefficient  matrices  [m^^] ,  [c^^]  and  [k^^] 
the  type  of  excitation  is  described  by  the  column  matrix 

In  summation  form.  Equation  1  may  be  written  as: 

^  Jdtjx^(t)  =  f„(t),  a  =  1,  2,  (2) 


*Multi  (or  single)  degree  of  freedom  systems  whose  motions  are  confined  along  a  single  dimensional  space  are  considered  one- 
dimcnsional.  Systems  which  arc  not  one-dimensional  are  consequently  both  multi-dimensional  and  multi-degree  of  freedom. 

is  the  damping  force  at  due  to  a  unit  velocity  at  is  the  elastic  force  at  x^^  due  to  a  unit  displacement  at  x^; 

contains  scalar  information  regarding  masses  as  well  as  first  and  second  moments  of  the  masses.  In  a  simple  system  consisting  of 
*  mass  particles  connected  by  springs  and  dash-pots,  having  no  coupling,  m^^^  would  be  either  the  mass  or  the  principal  moment  of 
inertia  of  the  ath  particle.  All  transfer  masses,  first  products  and  crossproducts  of  inertia  would  then  be  zero. 
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3.  MECHANICAL  IMPEDANCE  APPROACH 


Mechanical  impedance  is  formally  defined  for  linear  systems  over  the  domain  of  frequency  as 
the  ratio  of  the  Fourier  transforms  of  the  force  excitation  and  the  velocity  response.  By  Fourier 
transformation,  Equation  2  may  be  written  as  the  set 


N 

+k„^/ia>)  V^(a;)  =  (w)  ,  a  =  1,2,  ••'.N,  (3) 

/8=1 

where 

0)  -  2ni  =  circular  frequency, 

(oi)  =  Fourier  transform  of  velocity,  (t) , 

(cj)  =  Fourier  transform  of  force  excitation,  (t)  . 

Letting 

Za/9  (")  (4) 

represent  the  frequency-dependent  characteristic  property  of  the  system,  Equation  3  may  be 
written  as 


'a/3 


(5) 


or,  in  the  shorthand  matrix  form. 


M{v,}  =  K}'  («) 

In  these  equations,  are  complex  numbers  giving  the  ratio  of  the  transform  vibratory  force 
at  coordinate  a  to  the  transform  vibratory  velocity  at  coordinate  p.  They  are  functions  of  frequency 
and  may  be  called  impedance  parameters  of  the  system.  For  a  =  /3,  the  Z^^  are  called  point  im¬ 
pedance  parameters,  and  for  a  /  they  are  called  transfer  impedance  parameters.  Alternatively, 
Z^^  may  be  considered  as  the  ratio  of  vibratory  force  input  at  the  a  coordinate  to  vibratory  velocity 
response  at  the  /3  coordinate  when  all  other  coordinates  are  infinitely  restrained  (i.e.,  zero  veloci¬ 
ties).  Accordingly,  the  matrix  [z^^]  may  be  considered  as  the  mechanical  impedance  matrix  of  the 
system,  while  the  column  matrix  {v^}  represents  the  transform  of  velocities  corresponding  to  the 
transform  of  input  forces  } . 
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Heretofore,  the  developments  have  been  analogous  to  those  of  the  methods  of  dynamic  influence 
and  stiffness  coefficients.  Equation  6,  however,  expressed  in  terms  of  all  the  coordinates  which 
have  been  defined,  may  not  be  in  the  most  useful  form.  In  many  practical  situations,  some  of  these 
coordinates  are  not  of  interest  either  because  impedance  information  at  these  coordinates  is  not 
desired,  or  these  coordinates  are  not  accessible  for  physical  measurement.  Frequently,  the  co¬ 
ordinates  of  interest  are  those  at  which  subsequent  interconnection  with  another  system  is  antici¬ 
pated.  Consequently,  it  is  desirable  or  necessary  to  express  Equation  6  in  terms  of  matrix  quanti¬ 
ties  referred  to  only  the  coordinates  of  interest,  hereafter  called  "exterior  coordinates."  The 
coordinates  not  of  interest  are  called  "interior  coordinates." 

If  Equation  6  is  rearranged,  grouped  and  partitioned  according  to  "interior"  and  "exterior" 
applied  forces  such  that  n  coordinates  are  exterior  and  the  remaining  p  coordinates  are  interior. 
Equation  6  may  be  written  as 


[Zee]  [Zei] 

^{Ve}" 

"{FEr 

nn  np 

n  1 

^  ^ 

nl 

> 

[Zie]  [z„] 

{vj 

{Fx} 

_  pn  PP  _ 

^  pi  ^ 

The  subscripted  [z],  {v}  and  {f}  are  the  submatrices  corresponding  to  the  partition;  the  sub¬ 
scripts  E  and  I  denote  respectively  exterior  and  interior.  In  accordance  with  the  reciprocity 
theorem  for  linear  systems. 


[Zei]  =  [Zie]^ 


(8) 


where  T  denotes  transpose. 

The  expanded  form  of  Equation  7  provides  a  set  of  equations  from  which  we  can  develop  the 
matrices  entirely  in  terms  of  quantities  referred  to  exterior  coordinates.  Obtaining  the  submatrix 
{Vj}  from  the  second  equation  of  the  partitioned  form  and  substituting  in  the  first  equation  yields 


[Ze]K}  =  , 

nn  nl  nl 

where 

[^e]  ~  [^Ee]  “  [^El]  [^II]  ^  [^Ie] 

nn  nn  np  pp  pn 

{Fe'’}  =  [Zei]  [Zii]-*{Fi}  , 

nl  np  pp  pi 


and  [Zjj]  *  is  the  inverse  matrix  of  [Zjj] ,  which  is  assumed  nonsingular. 


(9) 


(10) 

(11) 
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Equation  9  may  be  termed  the  ’’exterior  equation”  of  the  system  since  all  quantities  now  refer 
to  exterior  coordinates.  The  matrix  [z^]  in  this  equation  then  represents  the  internal  impedance 
as  seen  from  the  exterior  points,  and  differs  conceptually  from  in  the  original  matrix  equa¬ 

tion  only  in  the  coordinates  to  which  it  is  referred.  Here  the  matrix  {f^}  is  defined  as  the  blocked 
force  matrix”  of  the  system,  and  it  represents  the  interior  excitations  as  transmitted  to  the  ex¬ 
terior  coordinates  for  exterior  coordinates  fixed  so  that  no  motion  occurs.  In  Equation  9,  the  neg¬ 
ative  sign  with  denotes  the  change  of  reactions  due  to  interior  forces  into  equivalent  externally 
applied  forces. 

If  a  system  contains  no  internal  source  of  energy,  the  matrix  of  interior  forces,  {^i}? 
consequently  the  matrix  of  blocked  forces,  {f^*^},  becomes  identically  zero.  Such  a  system  is 
termed  a  ’’passive”  system  as  opposed  to  an  ’’active”  system  when  {f^^}  is  nonzero.  According  to 
the  foregoing  developments,  the  activeness  or  passiveness  of  any  system  may  be  arbitrary  in  that 
it  depends  primarily  on  the  initial  definition  of  interior  and  exterior  coordinates. 

It  follows  from  Equation  9  that  the  matrix  of  exterior  velocities  {v^}  can  be  obtained  explicitly 
by  premultiplying  through  by  the  inverse  matrix  of  [z^.]  which  is  assumed  nonsingular, 

{V.}  =  •  (12) 

The  inverse  matrix  [Z]-'  is  termed  the  "matrix  of  mobility"  [Y] .  Frequently  the  use  of  [Y]  is  more 
desirable  than  [z]  (see  Section  7.3).  At  this  point  the  choice  is  trivial.  When  the  matrix  {Fj.}  of 
applied  forces  at  exterior  coordinates  is  zero,  the  velocity  matrix  of  Equation  12  is  significant 
enough  to  warrant  the  special  term  "free  velocity"  matrix,  expressed  by 

{Veo}  =  - 

In  terms  of  the  free  velocity  matrix.  Equation  12  becomes 

{Ve-Veo}  =  [Ze]-UFe}-  (13') 

An  alternative  form  of  Equation  9  may  be  obtained  in  terms  of  the  free  velocity  matrix  and  the 
mobility  matrix  (inverse  matrix  of  impedance)  of  the  system,  as  follows.  Premultiplying  Equation  6 
by  the  inverse  matrix  which  will  be  denoted  by  [Y^^]  yields 

K)  =  <*<) 


where  [Y]  is  the  mobility  matrix. 

Partitioning  in  terms  of  exterior  and  interior  coordinates,  Equation  14  can  be  written  as 


f 

[^ee]  [^ei] 

r{^En 

n  1 

nn  np 

nl 

< 

> 

{vj 

[^ie]  [^ii] 

{Fl} 

^  pi  > 

_  pn  PP  __ 

pi  > 
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Expanding  Equation  15  provides  a  set  of  equations  from  which  we  can  develop  the  matrices  entirely 
in  terms  of  quantities  referred  to  exterior  coordinates.  Obtaining  the  submatrix  {Fj}  from  the 
second  equation  of  the  partitioned  form  and  substituting  in  the  first  equation  yields: 

{V 

where 

[Ye]  " 

nn 

{^EO. 
nl 

The  determinant  of  [z^]  in  Equation  9,  or  of  [y^]  in  Equation  16  (when  damping  is  neglected) 
may  be  recognized  as  the  frequency  determinant  of  the  system  observed  from  the  exterior  coor¬ 
dinates.  Upon  expansion  it  yields  a  polynomial  in  the  roots  of  which  are  known  as  the  natural 
frequencies.  Cooresponding  to  each  value  of  an  amplitude  matrix  may  be  computed  for  each 
natural  mode. 

4.  MULTI-DIMENSIONAL  IMPEDANCE 

The  multi-dimensional  or  one-dimensional  aspect  of  a  system  is  characterized  by  the  nature  of 
the  matrices  of  Equation  7.  The  necessary  and  sufficient  condition  for  one-dimensional  systems  is 
uniformity  in  coordinate  types  (rectilinear  or  rotational)  which  define  excitation  and  response  of  the 
system  confined  to  a  single  dimensional  space.  For  example,  a  single  terminal  having  three  recti¬ 
lineal  and/or  rotational  types  of  coordinates  is  not  one-dimensional  within  the  scope  of  this  definition. 

The  general  nature  of  multi -dimensional  systems  invariably  requires  consideration  of  both 
rectilinear  and  rotational  types  of  impedance  analysis  and  measurement.  Rearranging  and  group¬ 
ing  with  respect  to  rectilinear  and  rotational  impedance,  as  well  as  to  interior  and  exterior  coor¬ 
dinates,  and  partitioning  accordingly.  Equation  6  may  be  expressed  as 


^Eo}  [^e]  {^E  }  ’ 

nl  nn  nl 

(16) 

[^Ee]  ”  [^El]  [^II]  ^  [^Ie]  ’ 
nn  np  pp  pn 

(17) 

=  [Yei]  [Yii]-*{Vi}  • 

(18) 

np  pp  pi 


[ZeeI 

[Ze?] 

[Ze'/] 

[Ze?] 

njn, 

"l"2 

"iPl 

"1P2 

np 

np 

[^ee] 

[Ze^I 

[z.1‘1 

[z//] 

{V} 

”2"  1 

n  2  ^  2 

"2PI 

"  2p  2 

n  2I 

^ 

n  2I 

[Zxf] 

[z/l] 

[Zif] 

[z//] 

{f/} 

Pin  1 

Pl"2 

PJPl 

P1P2 

Pjl 

pp 

[Z/e*] 

[z/eI 

[Z//] 

[Z//] 

{v/} 

__P2"l 

P2"2 

P2P1 

P2P2_ 

(19) 
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■where  the  superscripts  S  and  6  respectively  represent  rectilinear  and  rotational.  It  is  apparent  in 
Equation  19  that  there  are  four  possible  types  of  impedance: 


[z®*]  =  point  or  transfer  impedances  defined  by  ratios  of  rectilinear  forces  to  rectilinear 
velocities, 

[z*®]  =  transfer  impedances  defined  by  ratios  of  rectilinear  forces  to  rotational  velocities, 

[z^*]  =  transfer  impedances  defined  by  ratios  of  couples  to  rectilinear  velocities, 

[z^®]  =  point  or  transfer  impedances  defined  by  ratios  of  couples  to  rotational  velocities. 

In  accordance  -with  the  reciprocity  theorem  for  linear  systems,  these  impedances  and  their  matrices 
possess  the  folio-wing  convenient  transpose  properties: 


11 

N 

[z4?  . 

[z//]  = 

[z/iT 

11 

CO  ” 

[Zb^]  = 

[z//r 

[Zei"]  = 

[z/h*]"  . 

[Zb^/]  = 

[z/bI"  • 

(20) 

5.  INTERCONNECTION  OF  SYSTEMS 


In  Section  3  it  was  shown  how  the  performance  equation  of  a  system  can  be  expressed  in 
terms  of  coordinates  of  anticipated  interconnections,  heretofore  defined  as  exterior  coordi¬ 
nates.  Let  the  set  of  these  exterior  coordinates  corresponding  to  {Vj.}  be  common  to  two 
systems  which  are  otherwise  separated.  Block  diagrams  representing  the  two  systems  are 
shown  in  Figure  la.  The  systems  will  be  connected  as  shown  by  Figure  lb,  and  their  com¬ 
bined  behavior  will  be  examined  in  terms  of  interfacial  critical  (resonant  and  anti-resonant) 
frequencies  and  response. 


The  individual  performance  of  the  systems 
prior  to  connection  can  be  described  with  ref¬ 
erence  to  Equation  9  by  the  set  of  exterior 
equations: 

[z*]{v'}  = 

nn  nl  nl 

[z”]  {v”}  =  (f”  -  Fb”} 

mm  ml 


(21) 

(22) 


SYSTEM  I 

^‘1  '  '1^ 
Fjl,/ 

- ^0 - 

SYSTEM  H 

[d,] 

®  I 

L — ir>’° - 

["■>] 

a.  Systems  before  interconnection. 


b.  Composite  system. 


where  the  superscripts  I  and  ii  designate  a 
specific  system.  The  subscript  E  is  dropped 


Figure  1— Block  diagrams  of  (a)  subsystems 
and  (b)  composite  system. 
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refer  to  exterior  coordinates.  The  impedance  matrices  [z^]  and  [z^^]  are  roughly  the  counterparts 
of  the  "driving  point"  impedances,  looking  back  from  the  interface  commonly  associated  with  one¬ 
dimensional  interconnected  systems. 

Assuming  the  identification  of  coordinates  of  system  I  to  be  physically  and  mathematically 
compatible  with  those  of  system  1 1,  then  m  n  in  Equations  21  and  22.  When  the  systems  are  com¬ 
bined  to  form  a  composite  system,  Equations  21  and  22  must  satisfy  a  force  equilibrium  condition 

{R}  = 

and  a  velocity  compatibility  condition* 

{V} 


|fI  +fII  -(rb'  +  Fb")} 


=  {v^}  =  {v»}  . 


nl 


(23) 


(24) 


where  the  nonsuperscripted  quantities  refer  to  composite  quantities  at  the  interface. 

Satisfying  these  conditions  in  Equations  21  and  22,  and  subsequently  summing  up  these  equa¬ 
tions  yields 


[H]  {V}  ^  {R} 

nn  nl  nl 

where 

[H]  =  [Z^+Z”] 

nn  nn 

Equation  25  is  termed  the  "composite  systems  equation,"  which  is  the  equation  of  equilibrium  in 
terms  of  the  interconnection  velocity  {V}  for  the  complete  system  and  the  matrix  [H] ,  called  the 
"interconnection  matrix." 

Thus,  in  the  calculation  of  the  interconnection  velocity,  the  complete  system  has  been  regarded 
as  an  assembly  of  the  subsystems  subjected  to  the  equivalent  loading  {R}  of  Equation  23.  The  sum¬ 
mation  implies  addition  of  the  corresponding  interconnection  reactions  for  intercon¬ 

nections  blocked  (i.e.,  blocked  forces),  while  +F“}  is  the  loading  matrix  for  externally  applied 
forces  on  these  interconnections. 

In  the  case  of  free  vibrations. 


(25) 


(26) 


{pi}  =  {f^^}  =  {f^^}  =  =  0  , 


(27) 


*The  nature  of  the  mechanical  connections  at  the  interconnection  (i.e.,  boundary  conditions)  are  implicitly  defined  by  the  force  equili¬ 
brium  and  velocity  compatibility  equations. 
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and  consequently 


{R}  =  {0}  . 


(28) 


and  Equation  25  reduces  to 


[H]  {V}  =  {0}  , 

a  set  of  linear  homogeneous  equations  in  V.  The  roots  of  the  frequency  determinant 


(29) 


A(H)  =  0 

correspond  to  the  natural  frequencies  of  the  composite  system 


(30) 

observed  at  the  interconnection. 


6.  ANALYTICAL  APPROACH 

In  order  to  evaluate  the  matrices  of  the  foregoing  sections,  analytical  and  experimental  pro¬ 
cedures  are  recommended  in  the  sections  that  follow.  The  application  of  these  procedures  in  con¬ 
junction  with  automatic  data  processing  techniques  will  greatly  enhance  the  computational  schemes 
required  by  theory. 

The  basic  procedure  for  analytical  determination  of  the  desired  matrices  discussed  in  the 
foregoing  sections  may  be  summarized  as  follows: 

(i)  Choose  a  set  of  appropriate  grid  points  for  the  system  and  identify  the  required  general' 

ized  coordinates.  In  terms  of  these  coordinates,  obtain  the  matrix  equation  of  motions 
from  which  the  coefficient  matrices  [k^^]  and  [c^^]  are  determined. 

(ii)  Determine  the  matrix  of  applied  forces  (or  moments)  {F} . 

(iii)  By  elementary  matrix  operations,  determine  the  submatrices  [Zgg] ,  [Zji],  [ZjeJj  [^u]’ 

>  {Ve}»  {vJ,  {Fe}  and  {f,}  from  [z„^],  {F„}and{v^}. 

(iv)  Determine  [z^.]  and  {f^^}  from  Equations  10  and  11. 

(v)  Determine  the  response  {Vg}  from  Equation  12. 

(vi)  Setting  the  damping  coefficients  to  zero  in  the  frequency  determinant,  determine  the 
critical  (resonant  and  anti -resonant)  frequencies  and  the  corresponding  natural  modes. 

If  the  effect  of  connecting  a  second  system  to  the  first  system  is  desired,  the  following  addi¬ 
tional  steps  are  required: 


(vii)  Repeat  steps  (i)— (vi)  for  the  second  system. 

(viii)  Determine  the  interconnection  matrix  [H]  from  Equation  26. 
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(ix)  Determine  the  loading  matrix  {R}  from  Equation  23. 

(x)  Determine  the  response  at  the  interconnection  from  Equation  25. 

(xi)  Neglecting  damping,  determine  critical  (resonant  and  anti -resonant)  frequencies  of  com¬ 
posite  system  from  Equation  30  and  compute  corresponding  natural  modes. 


7.  MEASUREMENT  APPROACH 

Experimental  measurement  of  the  impedance  matrix  [z^.]  requires  the  acquisition  of  a  com¬ 
bination  of  experimental  data  which  will  supply  sufficient  information  for  evaluating  the  matrix 
elements  of  [Zj.].  The  following  schemes  for  experimentally  determining  matrix  [Zg]  are  based 
on  the  fact  that  if  a  performance  equation  (Equation  9  or  16)  is  to  truly  describe  the  dynamic  be¬ 
havior  of  the  system,  matrix  [Zg]  must  be  compatible  with  the  set  of  velocities  {Vg}  and  forces  {Pg} 
measured  on  the  system  at  the  specified  coordinates. 

7.1.  Determination  of  Impedance  Parameters 

7.1.1  Active  Systems 

In  the  experimental  determination  of  impedance  parameters,  the  measurement  of  free  velocity 
Vq  is  more  readily  accomplished  than  that  of  blocked  force  .  This  is  due  to  the  fact  that  ideal 
physical  restraint  conditions  required  for  block  force  measurement  are  seldom  realizable.  By 
substituting  Equation  13  into  Equation  12,  the  performance  equation  for  active  systems  may  be 
written  in  modified  form  as 


[Z]{v-Vo}  =  {F}  ,  (31) 

where  the  subscript  E  has  been  dropped.  Hereafter,  the  absence  of  E  implies  that  all  coordinates 
are  exterior.  For  a  system  of  n  coordinates,  the  number  of  impedance  parameters  will  be  n^. 
Theoretically,  only  n(n  +  l)/2  number  will  have  different  values,  since  Z^^  =  Z^^,  for  (a  /  /3)  .  It  is 
highly  desirable,  however,  to  evaluate  all  parameters  since  this  will  provide  a  means  of  valida¬ 
ting  the  assumption  of  linearity  of  the  system  by  use  of  the  symmetric  condition  Z^^  =  Z^^  (a  /  /3). 
Since  the  column  matrix  of  free  velocities  must  also  be  determined,  a  total  of  n  +  l  tests  are  re¬ 
quired  at  each  frequency  to  evaluate  all  the  impedance  matrix  elements  and  the  elements  of  free 
velocity  of  Equation  31. 

Test  1:  With  the  system  suspended  on  its  natural  restraints,  or,  in  the  case  of  a  free  system, 
on  a  low  frequency  suspension  such  that  all  exterior  coordinates  are  unattached  (zero  load  im¬ 
pedance),  operate  the  system  and  determine  the  phasor*  of  all  resulting  free  velocities.  Denote 
these  free  velocities  by  (cj) }  where  the  superscript  (1)  represents  test  number. 


*The  term  "phasor,”  as  used  here,  implies  magnitude  and  phase  of  a  steady  state  sinusoidal  quantity  where  phase  measurement  is  made 
with  respect  to  an  arbitrary  fixed  reference.  If  a  quantity  is  not  steady  state  sinusoidal,  Fourier  transform  of  the  quantity  is  required 
for  determining  the  phasors. 
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Test  n  +  i;  With  the  system  suspended  as  in  Test  1,  operate  the  system  and  apply  an  external 
force  at  one  of  the  exterior  coordinates  (e.g.,  the  a'"  coordinate)  and  determine  the  phasor  of  all 
resulting  velocities  at  all  exterior  coordinates  and  the  phasor  of  the  applied  force.  Denote  these 
velocity  phasors  and  this  applied  force  phasor  respectively  by  {v(2)  (w)}  and  {f(2)  (a>)}  ,  where  all 
elements  of  {f<2)}  except  the  a*  element  are  zeros.  Repeating  this  test  n  times,  each  time  ap¬ 
plying  a  force  at  successively  different  coordinates  and  determining  the  phasor  of  all  resulting 
velocities  and  the  applied  force,  provides  a  set  of  experimental  data  from  which  the  following  re¬ 
sulting  matrix  equations  may  be  solved  simultaneously  for  the  elements  of  [Z]  at  any  specified 
frequency,  say  : 


[z(-o)]  (-0)  (-0)}  "  K)}  •  (32) 

where  the  test  number  k  =  2,  3,  •  •  •  ,  n  +  1. 

7.1.2  Passive  Systems 

For  passive  systems,  the  matrix  of  free  velocities  {Vq}  is  identically  zero  and  Equation  31 
reduces  to 


[Z]  -{V}  =  {F}  • 

Accordingly,  tests  of  the  type  outlined  under  Test  (n  +  1)  for  active  systems  (except  now  there  are 
no  free  velocities)  provide  a  set  of  experimental  data  from  which  the  resulting  matrix  equations 
(similar  to  Equation  32)  which  follow  may  be  solved  simultaneously  for  the  elements  of  [Z]  of  the 
passive  system  at  any  specified  frequency 


[z(c.o)]  K)}  =  (^0)}  -  (34) 

where  k  -  1,  2,  •  •  •  ,  n.  Since  passive  and  active  systems  are  conceptually  related  by  the  absence 
or  presence  of  an  appropriately  determined  free  velocity  or  blocked  force  matrix,  these  determ¬ 
inations  based  on  Equation  34,  if  supplemented  by  an  independent  determination  of  free  velocities, 
may  also  be  employed  for  active  systems. 

7.2.  Determination  of  Mobiiity  Parameters 

The  solutions  of  Equation  32  or  34  when  n  >  3  will  invariably  require  the  use  of  digital  com¬ 
puters.  A  more  practical  procedure  which  considers  the  computation  of  elements  of  the  inverse  of 
the  impedance  matrix  (mobility  matrix)  can  be  made  available  as  follows;  Equation  32  in  its  inverse 

form, 


(35) 


{yx  K)  K)}  =  [z(^o)]'‘  K)}  ' 
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where  [Y]  is  the  mobility  matrix,  and  expanded,  can  be  written  as  follows  (assuming  force  excita¬ 
tion  at  a  single  coordinate  to  be  taken  in  numerical  order  with  the  subscript  of  the  coordinates): 


V/'’  M  -  Vo^.V  K)  =  (6;J  +  0  +  0  +  •••  +  0 

("o)  -  Vo'.y  ("o)  =  K)  («o)  +  0  +  0  +  •  ••  +  0 

("0 )  -  Vo^.n’  (“0 )  =  F/2)  +  0  +  0  +  •  •  •  +0 

Vr*’  ("o)  -  Vo'V("o)  =  0  +  0  +  •••  +  Y,„  K) 

V/"^«  -  VW  (a.J  =  0  +  0  +  •••  +  Y,„  (coJ 

=  0  +  0  +  •••  +  Y„„  {a,,)  .  (36) 

From  Equations  36  the  elements  of  the  mobility  matrix  are  readily  computed  as 

("^oj  (cOq)  ’ 


where  a,  yS  =  1,  2,  3,  •  •  •  ,  n. 

For  the  case  of  passive  systems  wherein  no  free  velocities  are  present,  Equation  36  and  con¬ 
sequently  Equation  37  reduces  to 


(38) 


where  the  superscripts  (/3  +  i)  have  been  modified  to  be  compatible  with  the  test  number  shown  in 
Equation  34. 

The  complete  determination  of  the  impedance  or  mobility  matrix  requires  performing  the  above 
procedures  over  all  frequencies  of  interest.  The  impedance  matrix  [Z]  may  be  obtained  by  per¬ 
forming  the  inverse  of  the  mobility  matrix,  i.e.,  [Z]  =  [Y]"^ 

7.3.  Significance  of  Mobility  and  Impedance  Parameters 

It  is  evident  in  the  foregoing  analytical  procedures  that  the  significance  of  mobility  parameters 
versus  impedance  parameters  is  primarily  one  of  utility.  Since  they  are  mathematically  related 
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by  the  matrix  inverse  relation 


[z]  =  [Y]"‘  , 

the  choice  may  depend  on  whether  force  or  velocity  response  is  considered  as  the  unknown,  or  may 
depend  on  which  is  the  more  desirable  parameter  for  describing  one’s  interest.  In  problems 
dealing  with  creating  motions,  the  use  of  mobility  may  be  preferred;  in  problems  dealing  with  im¬ 
peding  motions,  the  use  of  impedance  may  be  preferred. 


A  greater  significance  which  is  related  to 
the  measurement  of  these  parameters  is  dem¬ 
onstrated  by  the  invariant  and  variant  property 
of  the  elements  of  the  mobility  and  impedance 
matrix,  respectively  (Reference  4).  Consider 
a  system  of  three  degrees  of  freedom  as  shown 
in  Figure  2a.  The  performance  equation  of 
this  system  in  terms  of  the  mobility  matrix  is 


where  {V}  is  the  column  vector  of  resultant 
velocities  corresponding  to  the  column  vector 
of  applied  forces  {F}.  The  elements  of  the 
mobility  matrix  [Y]  may  be  evaluated  by  the 
method  of  Section  7.2.  When  only  is  applied 
and  Xj  is  considered  the  exterior  coordinate 
(Figure  2b),  the  performance  equation  is  simply 


a.  Block  representation  of  a  one -dimensional/  three  degrees 


of  freedom  system. 


b.  Block  representation  when  only  and  are  considered. 


ore  considered. 


Figure  2— Block  representation  of  a  system  with 
varying  number  of  coordinates  under  consideration. 


from  which 


=  Y  F 
^11  ^1 


The  performance  equation  when  only  F^  and  F2  are  applied,  and  and  X2  are  considered  the  ex¬ 
terior  coordinates,  (Figure  2c)  is 
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from  which  the  elements  of  [Y] ,  evaluated  by  method  of  Section  7.2,  are: 


Y 


11 


It  is  noted  that  Y^^  remains  the  same  as  when  only  is  considered.  Considering  all  three  co¬ 
ordinates  x^,  X2  and  X3  with  applied  forces  Fj,  F2  and  F3,  the  performance  equation  becomes 


where  the  elements  Y^^  (a,  /?  =  1,  2)  remain  the  same  as  when  only  x^  and  X2  are  considered  and 
the  additional  elements  are  evaluated  as 


Accordingly,  the  elements  of  the  mobility  matrix  for  the  system  are  seen  to  be  invariant,  in  this 
manner,  with  respect  to  the  number  of  coordinates  considered.  Now  consider  the  computation  of 
the  elements  of  the  impedance  matrix  using  Equation  39.  A  general  element  may  be  computed 
from 


■  A(Y) 


(41) 


where  the  numerator  and  denominator  are  respectively  the  adjoint  and  determinant  of  the  mobility 
matrix.  As  the  numerator  and  denominator  in  Equation  41  are  dependent  on  the  order  of  the  matrix 
[Y]  and  are  in  general  not  equal,  the  elements  of  the  impedance  matrix  consequently  are  variant 
with  respect  to  the  number  of  coordinates  considered.  This  implies  that  as  additional  coordinates 
are  considered,  each  previous  matrix  element  must  be  recomputed. 

The  invariant  property  of  the  mobility  matrix  makes  it  desirable  and  natural  to  measure  and 
use  mobility  parameters,  and  to  consider  impedance  parameters  as  derived  quantities  with  certain 
utilities.  Furthermore  the  use  of  the  term  "driving  point,”  frequently  associated  with  one¬ 
dimensional  systems,  possesses  true  counterparts  in  multi-dimensional  systems  when  used  only 
with  elements  of  the  mobility  matrix. 

For  active  systems,  the  choice  of  free  velocity  or  blocked  force  measurements  requires  some 
consideration.  Mathematically,  the  free  velocity  matrix  is  simply  related  to  the  blocked  force 
matrix  through  Equation  13.  Computationally,  free  velocities  can  be  obtained  from  blocked  forces 
if  the  mobility  matrix  is  known.  Experimentally  however,  the  measurement  of  blocked  forces  in 
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general  is  exceedingly  difficult  because  of  the  physical  constraint  required  by  theory.  The  meas¬ 
urement  of  free  velocities  is  more  easily  accomplished  and  consequently  the  more  natural  quan¬ 
tities  to  measure. 


8.  DETERMINATION  OF  IMPEDANCE  IN  PRESENCE 
OF  PARTIAL  EXTERNAL  RESTRAINTS 

In  the  laboratory,  the  impedance  matrix  [z^.]  of  a  system  sometimes  cannot  be  determined 
readily  for  ideal  free  conditions  (when  the  system  is  naturally  free)  or  for  natural  restraints* 

(when  the  system  is  naturally  restrained).  Some  of  the  free  coordinates  may  be  restrained  to  some 
degree,  or  some  of  the  natural  restraints  may  be  altered.  In  some  cases  the  required  measurement 
instrumentation  may  impose  significant  restraints.  A  formal  method  to  determine  the  impedance 
[Zj,]  of  the  free  system  from  an  impedance  [Zg]  determined  with  the  system  unnaturally  restrained 
in  a  known  fashion  is  developed  in  the  following. 


Consider  an  active  system  having  a  class  of  u  exterior  coordinates  which  are  free  of  restraints 
and  another  class  of  v  exterior  coordinates  which  are  restrained  in  a  known  fashion.  When  the  set 
of  free  coordinates  is  denoted  by  subscript  f  and  the  set  of  restrained  coordinates  by  r,  the  system 
when  entirely  free  of  restraints  (Figure  3a) 


may  be  expressed  as 


[Yff]  [Yf,] 


[Yh]  [Y.,] 


which  is  obtained  by  substituting 


[Z]"^  =  [Y] 

in  Equation  12,  and  partitioning  in  terms  of 
submatrices  corresponding  to  free  and  re¬ 
strained  coordinates.  When  the  restraints  are 
expressed  in  terms  of  a  known  impedance 
matrix  [zj  as  shown  in  Figure  3b,  the  applica¬ 
tion  of  the  composite  equation  (Equation  25)  at 
the  restrained  coordinates  yields 


Active  ^ 

System  j. 

^  -  - 

Active 

Restraint  System 
e 

2, 

Vr 

[zj 

b.  Block  representation  of  restrained  system. 

Active  ^ 

- 0 - ^  ,  Vf 

System  r 

— 0 — ►F,  + 

Vr 

F^  -  Z  V 

c.  Equivalent  block  representation  of  Figure  3b. 


Figure  3— Block  representation  of  a  system  under 
[a)  natural  restraints,  and  (b,  c)  unnatural  restraints. 


*Hereafter,  no  distinction  will  be  made  between  free  and  naturally  restrained  systems. 
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where 


{R}  - 

and  (v^  }are  the  interconnection  velocities.  As  shown  in  Figure  3b,  {f/  +F^'}  are  the  externally 
applied  forces  at  the  restraints. 

The  experimental  evaluation  of  [Y]  of  the  restrained  system  which  will  be  denoted  by  [Y]  may 
be  accomplished  in  a  manner  similar  to  that  of  Section  7.2  by  applying  components  of  {f^}  and 
{f^'  -f-Fj}  one  component  at  a  time,  measuring  the  resulting  velocity  at  all  coordinates,  and  finally 
solving  for  the  elements  of  [Y]  in 


{v,r 

Iv„]  [V,] 

ul 

>  = 

UU  uv 

ul 

> 

U 

l> 

[v„l  P„1 

{R} 

.  J 

I 

vu  vv 

(45) 


The  objective  is  to  determine  the  unknown  [Y]  of  the  "free"  system  from  the  "measured"  [Y]  of 
the  restrained  system  and  the  known  restraint  impedance  [z^] .  This  may  be  accomplished  by  con- 
sidering  the  effect  of  the  restraint  system  as  reaction  forces  on  the  free  system  as  shown  in  Fig¬ 
ure  3c  when  {F/}  and  (f^'}  are  respectively  the  externally  applied  forces  on  the  free  and  the  re¬ 
straint  system  prior  to  connection,  and  {f^}  and  [z^]  {v^}  are  the  reaction  forces  due  to  the 
restraint  system.  Accordingly,  the  net  forces  acting  at  the  r  coordinates  of  the  free  system  of 
Figure  3c  are 


{fJ  =  {f/+f;}  -  {F^'’}  -  [zj{vj  • 


Using  Equation  44, 


{F^}  {R+F^^}  -  [ZJ{VJ  • 


{R}  =  {F,  -F,^}  +  [ZJ{VJ  • 

The  substitution  of  Equation  48  in  Equation  45  yields 


[Y„]  [Yj 

{Ff-r,*’} 

ul  1 

1  .a 

UU  U  V 

1> 

[V,,]  [?„) 

{F^-F,»}  +  [ZJ{VJ 

L  VI  J 

vu  vv 

(46) 


(47) 


(48) 


(49) 
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The  determination  of  [Y]  from  [Y]  may  be  accomplished  by  suitably  comparing  Equation  49 
with  Equation  42.  This  is  successfully  performed  by  first  noting  that,  as  the  result  of  the  equiva¬ 
lence  of  Figures  3b  and  3c, 


{V}  =  {V}  ’ 

and  then  operating  on  Equation  49  to  yield  the  form  of  Equation  42.  This  operation  yields  the 
desired  relation: 


[Yff]  [Yfr] 


[I] 


uu 


[[Vj  -  [Y„]]-' 


[Vr,]  [V,] 


[V.] 


[0] 


vu 


[».]  [K1  -  [v„l]“ 


[Y.f]  [y„] 


(50) 


Frequently  it  is  desirable  to  determine  the  effects  on  the  system’s  mobility  matrix  as  a  result 
of  some  known  partial  restraints.  This  is  accomplished  by  solving  [Y]  in  terms  of  [Y]  in  a  manner 
similar  to  the  above.  Accordingly  the  result  is: 


[Y„]  [Y,J 

UU  uv 

[I] 

uu 

-  [Y,J  [[yJ  +  [y„]]'‘ 

U  V  vv 

[Yff]  [YfJ 

UU  UV 

(V„l  (v„] 

vu  vv  _ 

[0] 

vu 

[Yj  [[yJ  +  KJ]" 

v  v  vv  _ 

[Y.r]  [Y,J 

VU  vv  _ 

If  all  coordinates  are  restrained.  Equation  51  reduces  to 


(51) 


=  KJ-'  +  [Y.]-'  . 

a  form  of  the  Interconnection  Matrix  previously  defined  by  Equation  26.  The  matrix  of  impedance 
[Z]  or  [Z]  may  be  obtained  by  application  of  the  matrix  inverse  relationship  between  the  imped¬ 
ance  and  mobility  matrix,  i.e. 

iZff]  [Zfr] 

Kr]  [Z.J 


[^ff]  [^fr] 

[z.f]  [z„] 


[Yff]  [Y,J 

[^rf]  [Y„] 


[I] 


[y,.] 

[Y..]  [Y„] 


=  [I] 
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9.  DETERMINATION  OF  FREE  VELOCITY 

IN  PRESENCE  OF  EXTERNAL  RESTRAINTS 


Frequently  the  free  velocity  matrix  {v^q}  of  an  active  system  cannot  be  determined  readily 
under  ideal  free  conditions.  It  is  therefore  desirable  to  develop  a  formal  relationship  to  determine 
the  free  velocity  matrix  of  the  system  from  the  velocity  matrix  determined  with  the  system  re¬ 
strained  in  a  known  fashion. 


As  in  the  previous  section,  the  desired  relationship  may  be  obtained  by  considering  some  of 
the  exterior  coordinates  f  to  be  free  and  the  remaining  coordinates  r  to  be  restrained  in  a  known 
fashion  (Figure  3).  Thus,  when  the  system  is  free.  Equation  16  describing  the  free  velocity  of  the 
active  system  may  be  written  as 


-  {Vo.j" 

U  1  U  1  1 

[Va] 

uu 

[v.r 

uv 

ul 

^  vl  vl  ^ 

>  - 

K,] 

vu 

[Y.], 

vv 

{F.}' 

vl  J 

which  is  obtained  by  partitioning  in  terms  of  submatrices  corresponding  to  free  and  restrained  co¬ 
ordinates.  When  the  system  is  restrained  by  impedance  [z^]  as  shown  in  Figure  3b,  the  effect  of 
the  restraint  system  on  the  free  system  may  be  considered  in  terms  of  reaction  forces  as  shown 
in  Figure  3c.  Thus,  the  net  forces  acting  on  the  free  system  at  the  restrained  coordinates  are 
given  by  Equation  46.  Accordingly,  replacing  {f^}  in  Equation  52  by  Equation  46,  and  denoting  the 
measurement  of  velocities  {v^}  and  {v^}  under  the  restrained  conditions  as  {v^  }  and  {v^}  respec¬ 
tively,  yields 


(V,}  -  {».,,) J 

For  free  velocity  measurements  at  the  free  and  restrained  coordinates, 

{Ff}  =  {F/}  =  {f;}  =  0  . 

Substituting  the  condition  of  Equation  54  in  Equation  53  and  simplifying  yields 
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(53) 


(54) 


(55) 


which  is  the  desired  relation  of  free  velocities  in  terms  of  velocities  measured  under  the  restraint 
condition.  The  procedure  consists  of  first  determining  the  mobility  of  the  free  system  [Y]  by  the 
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method  of  Section  8  and  the  blocked  forces  of  the  active  restraint  system.*  Next,  measure  the 
velocities  of  the  active  system  at  the  free  and  restrained  coordinates  (v,}  and  {vj  respectively 
and  apply  Equation  55  to  obtain  the  true  free  velocities  {V(, ,}  and{V(|  j  of  the  active  system  using 
the  known  mobility  [Y] .  It  may  be  seen  that  if  all  coordinates  are  free  ([yJ  ^ 


and 


as  would  be  expected. 


10.  EXAMPLES  OF  APPLICATION 
10.1.  Analytical  Approach 

Consider  the  two  dimensional  system  of  Figure  4  in  which  two  bodies  Mj  and  M2  are  supported 


by  elastic  and  damping  elements  located  at  their 
cal  with  respect  to  the  axis  passing  through  the 
center  of  gravity  of  the  bodies.  Vibration  of  the 
system  results  from  an  oscillating  force  act¬ 
ing  upon  Mj.  The  force  is  applied  at  a  point  on 
the  plane  of  symmetry,  at  a  distance  e  above 
the  center  of  gravity  ofM^  and  has  a  horizon¬ 
tal  component  coscot,  and  a  vertical  component 
f  y  sin  cot  . 

The  degrees-of-freedom  of  the  center  of 
gravity  are  as  shown  in  the  upper  right  corner 
of  each  body,  and  for  small  distortions  about  the 
system's  stable  equilibrium,  these  are  taken  as 
generalized  displacements  (/3  =  1,  2,  •  •  • ,  6). 
In  this  example  the  coordinates  x^,  x^,  and  x^ 
will  be  referred  to  as  "exterior." 

10.1.1  Energy  Expressions 

The  kinetic,  potential  and  dissipation 
energy  of  the  system  are  given,  respectively. 


corners,  as  illustrated.  The  system  is  symmetri- 


Figure  4— Lumped  parameter  equivalent  of  a 
simplified  active  two-dimensional  system. 


*Note  that  {f^}  can  be  obtained  from  free  velocity  measurements 


using  Equation  13* 
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=  M,  (i,2 +^2^^  +  Ij  +  Mj  (i^2 +i/)  +  Ij  , 

2Kj  Xj^  +  2K^  b^  x^^  +  2K2  X2^  +  2K2  a^  x^^ 

+  2K3  (Xj  -  x^Y  +  2K3  (CX3  -dXg)2 

+  2K^  (Xj-Xj)^  +  2K^  (0X3 -exg)^  , 

=  2Cj  Xj^  +  2Cj  X3^  +  2C2  +  2C2  x^^ 

+  2C3  (^1-^4)^  +  2C3  (0^3 

+  2C4  +  2C4  (eij-eig)^  •  (56) 

10,1,2  Equations  of  Motion 

Applying  Lagrange equation  for  nonconservative  systems  yields  the  equations  of  motion: 

Ml  Xj  +  2(Ci  +  C3)  Xj  -  2C3  x^  +  2(Ki  +  K3)  Xj  “  2K3  x^  -  f  1  , 

Mj  X2  +  2(C2+C^)x2  “  20^X3  +  2(K2+K^)x2  -  2K^  Xg  =  £2  , 

Ij  ^3  +  2(Cjb2  +C3  +C3C2  +C^  62)^3  -  2(03  cd  +  C^e2) 

+  2(Ki  b^+K2a^+K3C^+K^e^)x3  -  2(K3  cd  +  K^e^)xg  = 

^2^4  '203X1  +  2C3  X4  -  2K3X1  +  2K3X4  =  £4  , 

M2  Xg  -  2C4X2  +  2C4  ^5  -  2K4X2  +  2K4  X3  =  fg  , 

^2  ^6  '  2(03  cd  +C4  e^]  X3  +  2(03  d^  +C4  Xg 

-  2(K3  cd +K4  e2)  X3  +  2(K3  d2 +K4  e2)  Xg  =  fg»  (57) 

where  £^  (a  =  1,  2,  •  •  •  ,  6)  are  generalized  forces  to  be  determined  later. 


■ 


2T 


2P 


2D 
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10.1.3  Coefficient  Matrices 

From  the  equations  of  motion,  the  coefficient  matrices  are  obtained  as: 


mil  *^22  ^1 


^33  ’ 


*^44  *^55  ^2  ’ 


‘66  2 


=  0  ;  a  /  /3  =  1,2,  *  *  ■  ,  6) 


ku  =  2(K,  +K; 


k,,  =  2{K,+Kj  , 


k,,  -  2(Kjb2+K,a2  +K3c2+k^e2) 


k44  =  2K3 


kjs  =  2K,  , 


=  2(K3d^+K,e  =  ) 


^14  ”*  ^41 


^25  ^52  ^4 


^36  =  k„  =  -2(K3cd  +  K,e' 


All  others  are  zero. 
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Cji  2(Cj  +C3)  , 

^22  “  ^(^2  '''^4)  ’ 

C33  "  2(c,b2+C,a^+C3c2+C,e2) 


■44  ^^3 


^55  ^^4 


C66  =  2(C3d^+C,e^)  . 

^14  "  ^41  “  “  ’ 


^25  *^52 


2(C3  cd  +  C.e^) 


All  others  are  zero. 

From  Equation  4, 

Zn  =  icoMj  +  2(C, +C3)  +  2(K, +K3)/ia,  , 

Z33  =  icoM,  +  2(C3  +  cJ  +  2(K3  +  Kj/i«  , 

Z33  =  ia;Ij  +  2(C,b2+C3a2  +  C3c2+C^e2)  +  2(Kjb2+K3a2+K3c2+K3e2)/i<i) 


Z^^  =  itaMj  +  X3  +  2K3/icJ  , 


iwMj  +  2C^  +  2K^/ia)  , 


io)  Ij  +  2(03  d^  +C^  6=  )  +  2(K3  d^  +K4  e^yio)  , 


Z4.  =  -  2(C3+K3/i^)  , 


Z52  =  -2(C,  +Vic.). 


Zga  =  -  2(C3  cd  +  C^e^)  -  2(K3  cd  +  K^e2)/i^ 
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The  type  of  impedance  has  been  designated  on  the  extreme  right.  Accordingly,  rearranging  rows 
and  columns,  and  subsequent  partitioning  of  [z^^]  yields  the  following  submatrices:* 


10.1  A  Generalized  Forces 

The  total  virtual  work  done  by  the  applied  forces  on  the  system  is 

8W^  =  fy  8x2  +  8x1  ~  6  8x3  .  (63) 

♦These  submatrices  are  required  by  Equation  10  to  determine  [z^j- 
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The  generalized  force  (corresponding  to  is  defined  by 


(64) 


from  which  fj  =  =  f  ,  fj  -  -f^e.  Thus 


{F:} 


F 

X 

€F 


{Fe}  =  {0} 


(65) 


where  F^,  F^  and  -  eF^  are  the  Fourier  transforms  of  fj ,  and  £3  respectively. 

10,1.5  Determination  of 

In  Equation  10,  the  inverse  matrix  of  [z^j]  is  required  for  the  determination  of  [Zg]-*  ^et  us 
denote  the  elements  of  by  Then  is  related  to  by 


ak  ^k/3  ^a/3 


(66) 


where  is  the  so-called  ''Kronecker  delta.’' 

Since  [Zjj]  is  a  diagonal  matrix.  Equation  66  reduces  to 


Thus 


^Il.aa  baa 


1  . 


0 
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0 
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10.1.6  Determination  of  [z^  ] 

From  Equation  10,  the  general  element  of  [z^,]  is  given  by 


"E.ayS 


"EE,  a/3 


Lz: 

'£=1  k=i 


"El.ak  bii, 


(a,  /3  =  1,  2.  3) 


(67) 


(68) 


♦The  inversion  for  this  particular  example  is  easily  obtained  from  Equation  66.  In  general,  inversion  of  higher  orders  may  require 
machine  computation. 
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Since,  for  this  example: 


then 


and 


Thus 


^EE.ayS 

7 

^El,ak 


(o.  f  J3)  , 

(a  /  k)  , 

(k  /  -t)  , 
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10,1.7  Determination  of  {Fe%} 
From  Equation  11, 


3  3 


E,a 


k=l 


(a  =  1,  2,  3) 


but  for  this  example,  -  z^j  ^  ,  from  which 
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(71) 
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Figure  5— Lumped  parameter  equivalent  of  the  composite 
system  composed  of  systems  I  and  II,  with  connection  at 
center  of  gravity  of  . 


Figure  6— Lumped  parameter  equivalent  of  a 
simplified  passive  two-dimensional  system. 


10.1,8  Determination  of  {Vj,} 
From  Equation  12, 


r  ^ 

V4 

K,a}  = 


f-^,4  (Zn 

~  ^25  (^22  ^55  ^  ^25)  ^  ^ 

^6  (Z33 


(72) 


10.1.9  Effect  of  Subsystems 

Next  we  shall  consider  the  effect  of  connecting  a  second  system  to  the  exterior  coordinates  x^, 
Xg,  and  Xg  (Figure  5).  The  system  is  assumed  to  be  passive  (Figure  6). 

Letting  superscripts  I  and  II  denote  quantities  corresponding  respectively  to  the  first  and 
second  system,  the  following  is  obtained: 


where 
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+  gMc."  +Kn/ico)  . 

Z  »  =  C  n  .  K  , 

Zs"  =  cn  +  K^/i..  , 

Z»  =  f^(C."+K»/i.)  , 

Z  »  =  -  (C.n  .K -/ic.)  , 

Z3V  =  -  (C»  .K»/i.)  , 

ZaV  =  -  gf  (C”  +K»/i.,)  .  (73) 

From  Equation  26,  the  interconnection  matrix  is 


II 


-(za)’Aa] 


L  =  I 


II 


L=I 


(74) 


The  roots  of  A(H)  =  0  are  the  critical  frequencies  of  the  composite  system. 
If  0^^  denotes  the  elements  of  then 


[Qa,] 


l/H^^  0  0 

0  1/Hj2  0 

0  0  I/H33 


(75) 
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It  follows  from  Equation  23  and  25  that  the  response  of  the  composite  system  at  the  interconnection 
of  the  two  systems  is: 

V4  =  -  (ZH/Zn)F,/j[z44  ’ 

V5  =  -  (Z2sA22)f7£[Zs"-  (225^/^32]  ' 

/  L=I 

Ve  =  (Z3V233 )  «  Fx/  fy-l,  -  )  y Z3''3]  •  (76) 

/  L=I 

10.2.  Measurement  Approach 

The  primary  purpose  here  will  be  to  demonstrate  the  experimental  evaluation  of  the  mobility 
matrix  [y^  ]  of  each  of  the  systems,  to  obtain  their  corresponding  interconnection  matrix  and  to 
describe  the  resulting  response  at  the  interconnection. 

10.2,1  Determination  of  Mobility  Parameters 

Equation  35  for  passive  systems*  is 


{vcx)  (..„)}  =  [yK)]{F('<)  . 

where  k  =  3  since  there  are  three  exterior  coordinates  in  each  system. 


(35) 


10.2.2  System  I 

With  the  system  suspended  on  a  low-frequency  suspension  such  that  all  exterior  coordinates 
(x4,  X5,  xg)  are  unattached  (zero  load  impedance  or  infinite  mobility),  apply  a  sinusoidal  force  of 
constant  frequency  co^  at  x^  and  measure  the  phasor  (referencing  phase  angle  with  respect  to  the 
applied  sinusoidal  force)  of  the  resulting  velocity  at  x^  ,  x^  and  Xg .  Denote  these  velocity  phasors 
and  this  single  applied  force  phasor  respectively  by 
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(77) 


♦Although  system  I  is  truly  active,  the  use  of  Equation  35  for  passive  systems  supplemented  by  an  independent  determination  of  free 
velocities  is  equally  applicable. 
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Repeating  the  above  test  two  more  times,  each  time  applying  the  sinusoidal  force  of  frequency 
successively  at  and  x^,  provides  the  additional  set  of  measurements  similar  to  (77): 
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In  Equation  80  the  condition  (a  /  /3)  provides  a  check  on  the  linearity  of  the  system. 

The  complete  determination  of  the  mobility  matrix  requires  repeating  the  above  procedures 
for  all  frequencies  of  interest. 

10.2.3  System  U 

The  determination  of  the  mobility  matrix  for  system  II  is  accomplished  in  a  manner  similar 
to  that  for  system  I.  The  exterior  coordinates  of  this  system  are  defined  as  and  Xg^^  as 

shown  in  Figure  6. 

10.2.4  Determination  of  Free  Velocities 

In  this  example,  it  is  necessary  only  to  determine  the  free  velocities  of  system  I  since  system 
II  is  passive.  With  the  system  suspended  as  before,  the  resulting  free  velocity  phasor  at  all  ex¬ 
terior  coordinates  due  to  the  excitations  and  at  frequency  may  be  measured  directly.  A 
convenient  reference  for  phase  angle  measurements  is  the  free  velocity  at  one  of  the  coordinates. 
Assuming  that  the  phase  angle  of  these  measurements  is  made  relative  to  the  free  velocity  at  co¬ 
ordinate  x^,  the  free  velocity  phasors  may  be  expressed  as 
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(81) 


The  complete  determination  of  free  velocity  requires  operating  the  system  over  all  frequencies  of 
interest  and  performing  the  above  procedure. 
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10.2.5  Interconnection  Matrix 


If  we  denote  the  [y^]  for  the  active  system  and  passive  system  respectively  by  superscripts  I 
and  II,  the  interconnection  matrix  may  be  obtained  from  Equation  26  as 

[h(c<,)]  =  [yi  (ao]-i  +  [y"  (co)]-»  . 

where  the  matrix  inversion  of  [Y(aj)]  may  require  digital  machine  computations. 

10.2.6  Prediction  of  Composite  Response 

The  response  at  the  interconnection  of  the  two  systems  may  be  predicted  by 

{V(c.)}  =  -  [H(<.)]->{f^\co)}  .  (83) 

where  {f'’’}  can  be  obtained  from 

=  [y'  (")]■' K  (")}  •  (84) 


11.  DISCUSSION 

The  objective  of  this  paper  is  to  present  a  systematic  approach  for  the  analytic  and  experi-^ 
mental  determination  of  mechanical  impedance  and  its  associated  effects  on  dynamic  response  of 
complex  mechanical  systems.  The  nature  of  the  formulation  of  the  approach  is  sufficiently  general 
for  extension  to  cases  of  many  types  of  distributed  systems  subject  to  varied  steady  state,  transient 
and  random  excitations. 

Although  the  concepts  related  to  the  use  of  mechanical  impedance  expressions  have  been  ex¬ 
plored  extensively  in  the  development  of  methods  for  shock  and  vibration  control,  they  have  been 
limited  to  a  large  extent  to  systems  which  permit  analysis  of  only  one -dimensional  types  (Refer¬ 
ence  2).  Limited  and  specialized  application  to  multi -dimensional  systems  has  been  studied  only 
recently.  As  shown  in  this  study,  the  treatment  of  mechanical  impedance  of  multi -degree  of  freedom/ 
multi-dimensional  systems  may  be  approached  by  adopting  certain  aspects  of  matrix  methods  of 
structural  analysis,  together  with  appropriate  interpretation  of  quantities  related  to  the  concepts 
of  mechanical  impedance  (e.g.,  blocked  force,  free  velocity,  point  and  transfer  impedance). 

The  problem  of  describing  complex  mechanical  systems  so  as  to  consider  the  information 
available  about  the  individual  subsystems,  involves  the  formulation  of  impedance  matrices  and 
interconnection  equations  wherein  the  conditions  of  force  equilibrium  and  velocity  compatibility, 
implicitly  imposed,  are  equivalent  to  physical  connection  of  the  subsystems.  The  characteristic 
values  of  the  interconnection  matrix  correspond  to  critical  frequencies  of  the  composite  systems  as 
observed  at  the  interconnection.  Moreover,  the  methods  of  interconnection  developed  are  applicable 
to  many  systems  of  complex  combinations.  Systems  of  mixed  types  such  as  electro- mechanical 
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and  acoustic-mechanical  may  be  treated  using  these  methods.  With  regard  to  environmental  test¬ 
ing,  the  methods  of  interconnection  further  provide  a  sound  theoretical  basis  for  performing  sub- 
assembly  testing  which  would  obviate  the  demand  for  increasingly  large  complex  environmental 
test  facilities  for  future  generation  aerospace  systems. 

On  the  basis  of  the  results  of  this  study,  the  following  recommendations  are  made: 

1.  Suitable  experiments  should  be  undertaken  to  determine  the  proper  extent  of  the  practical 
application  of  the  methods  herein  presented. 

2.  Effort  should  be  directed  toward  developing  a  computer  program  for  mathematical  compu¬ 
tation  of  the  required  matrix  relationships. 

3.  Effort  should  be  directed  toward  the  development  of  sensors  and/or  data  reduction  methods 
for  the  measurement  of  multi-dimensional  impedance  parameters. 

4.  Effort  should  be  directed  toward  the  investigation  of  experimental  techniques  for  determin¬ 
ing  these  parameters  using  random  and  transient  types  of  force  excitations. 

5.  The  feasibility  of  applying  the  results  of  this  study  to  environmental  test  simulation  should 
be  investigated  in  further  studies.  Included  should  be  a  laboratory  study  of  simulation  and 
control  of  impedance  of  multi-dimensional  models. 
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Appendix 

Symbols 


a,  b,  c,  d,  e,  f,  g 
C 
D 

e 

h 
SS 
5(9,  eh 

E 

EE 
El,  IE 

fa 

fx 

fy 

I 

II 

I 

K 

M 

N 

P 

Fa(^) 


linear  dimensions  for  example  problem 
damping  coefficient 
dissipation  energy 
distance  above  c.g.  of  mass 
superscript  for  rectilinear  quantities 
superscript  for  rectilinear  couplings 

superscripts  for  rectilinear— rotational  couplings,  or  conversely 
subscript  for  exterior  quantities 
subscript  for  exterior  couplings 

subscripts  for  exterior-interior  couplings,  or  conversely 

generalized  force  corresponding  to  ath  generalized  coordinate 

applied  force  in  x -direction 

applied  force  in  y-direction 

subscript  for  interior  quantities 

subscript  for  interior  couplings 

mass  moment  of  inertia 

stiffness  coefficient 

mass 

number  of  degrees-of-freedom 

potential  energy 

Fourier  transform  of  (t) 
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(oj)  Fourier  transform  of  (t) 
n  number  of  exterior  coordinates 
0  subscript  for  rotational  quantities 
66  subscript  for  rotational  couplings 
T  kinetic  energy 
t  time  variable 

I,  II  superscripts  for  specific  subsystems 
w  circular  frequency 

displacement  corresponding  to  /3th  coordinate 
x^  velocity 
x^  acceleration 

(6;)  frequency  dependent  characteristic  property  of  system 
8W^  total  virtual  work  done  by  applied  force 
8x^  virtual  displacement 
{  }  denotes  column  matrix 
{  denotes  transpose  of  {  } 

[  ]  denotes  square  or  rectangular  matrix 
[  ]  ^  denotes  transpose  of  [  ] 

[  ]"^  denotes  inverse  matrix  of  square  matrix 
inertia  matrix 
damping  matrix 
stiffness  matrix 
I transform  force  matrix 
{f^^}  blocked  force  matrix  (transform) 

{v^}  transform  velocity  matrix 
{Veo}  transform  free  velocity  matrix 
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[z^^]  mechanical  impedance  matrix  for  system 
[Zj.]  impedance  matrix  looking  back  from  exterior  coordinates 
[Y]  inverse  matrix  of  [Z] 

[zj  impedance  matrix  of  restraint  system 
[Ha^s]  interconnection  matrix 
A(H)  determinant  of  [H^^] 

[q^^]  inverse  matrix  of  [H^^] 
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*'Th€  aeronautical  and  space  activities  of  the  United  States  shall  be 
conducted  so  as  to  contribute  ...  to  the  expansion  of  human  knowl¬ 
edge  of  phenomena  in  the  atmosphere  and  space.  The  Administration 
shall  provide  for  the  widest  practicable  and  appropriate  dissemination 
of  information  concerning  its  activities  and  the  results  thereof.*^ 

—National  Aeronautics  and  Space  Act  of  1958 
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